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Abstract 

This work treats the problem of weighted least squares fitting of a 3D 
Euclidean-coordinate transformation matrix to a set of unit vectors measured m the 
reference and transformed coordinates. A closed-form analytic solution to the problem 
is re-derived. The fact that the solution is the closest orthogonal matrix to some 
matrix defined on the measured vectors and their weights is clearly demonstrated. 
Several known algorithms for computing the analytic closed form solution are 
considered. An algorithm is discussed which is based on the polar decomposition of 
matrices into the closest unitary matrix to the decomposed matrix and a Hermitian 
matrix. A somewhat longer improved algorithm is suggested too. A comparison of 
several algorithms is carried out using simulated data as well as real data fro™ .the 
Upper Atmosphere Research Satellite. The comparison is based on accuracy and time 
consumption. It is concluded that the algorithms based on polar decomposition yield a 
simple although somewhat less accurate solution. The precision of the latter 
algorithms increase with the number of the measured vectors and with the accuracy of 
their measurement. 
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1. INTRODUCTION 


The problem of attitude determination from vector observations is as follows. A 
sequence, bj, i=0,l,2,...K of unit vectors is given. These unit vectors are the 

result of measurements performed in vehicle body axes of the directions to known 
objects. The sequence, r, i=0,l,2,...K of unit vectors, is the sequence of the 

corresponding representation of these directions with respect to some reference 
coordinate system. We wish to find the attitude matrix, A, such that the cost 
functional p(A) defined as follows 

K 

P (A) = Z 2 aj I I bj - Arj I 1 2 (1) 

i = 1 

is minimized. This problem, which is basically a least-squares fit problem for the 
attitude matrix, A, was posed in [1] and is generally known as Wahba’s problem. This 
problem has been treated extensively [see, e.g. 2-11]. 

In the next section we derive an analytic solution to Wahba’s problem, then in 
Section III we show, in a rather simple way, that this solution is actually the 
closest orthogonal matrix to a matrix defined on the reference and measured unit 
vectors and bj respectively, and on their relative weight. Several algorithms for 

computing the attitude matrix are considered in that section. The connection between 
polar decomposition of matrices and the solution to Wahba’s problem is then discussed 
in Section IV. Two algorithms for computing the solution, which are based on the 
polar decomposition, are considered. A numerical comparison between these algorithms 
and other suggested ones, using simulated as well as real satellite data, is 
presented in Section V. The conclusions of this work are finally presented in Section 


II. DIRECT SOLUTION OF WAHBA’S PROBLEM 

Since only the relative value of the weights, a^, matter, we may, with no loss 
of generality, normalize the weights to give 

K 

E a i = i 

i = 1 

It can be shown [2] that 


p(A) = 1 - tr( AB 1 ) 


where tr denotes the trace of a matrix and 

K 


B = V a-b-r. 
L ill 

i = l 


( 2 ) 


(3) 


We seek the orthogonal matrix, A, which minimizes p(A). Obviously, that matrix 

maximizes tr(AB ). Using the method of Lagrange multipliers, we can incorporate the 

orthogonality constraint on A into the maximization problem of tr(AB T ) by defining 
the new functional p*(A) 
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(4) 


p*(A) = tr(BA T ) + trtjLCAA 7 - I)] 

where I is the 3x3 identity matrix. The matrix L is a matrix of Lagrange multipliers 
scaled to enable the inclusion of the one half factor which is added for simplicity 
of the ensuing derivation. Also note that with no loss of generality, we may choose L 
to be symmetric. The new cost function, p*(A), can be written as follows 

p*(A) = tr[(BA T ) + j L(AA T - I)] (5) 

Use now the directional derivative to maximize p*(A). To accomplish this, express A 
as follows 

A = A + eH (6) 

o 

where A is the A matrix which maximizes p*(A), e is a scalar variable, and H is any 
o . . . 

3x3 real matrix. Note that A in (6) is expressed as a sum of the maximizing matrix, 

A and a "step", e, in the "direction" of H. Also note that any real 3x3 matrix can 
°’ 

be expressed in this way. Substitution of (6) into (5) gives 

p’(e) = tr{B(A Q + eH) T + j L[(A q + eH)(A Q + eH) T - I]} (7) 

Next differentiate p’(e) with respect to e to obtain 

d P d M = tr[BH T + \ LU(A J o + eH T ) + \ L(A q + eH)H T ] (8) 

A necessary condition for p’( e ) 1° have a maximum at A^ is 

d P'( e > = 0 for all H (9) 

^ e=0 

Applying (9) to (8) yields 

tr[(B + LA q )H T ] = 0 for all H (10) 

The latter can exist only if 


B + LA =0 
o 

or, assuming L is non-singular, 

A = l'b (ID 

o 

Using (11) in the orthogonality constraint on A q and making use of the fact that L is 
symmetric we obtain 

A A T = L-'bbTl 1 = I 
o o 

which yields 
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T 2 
BB 1 = L z 


T 

The matrix BB is a positive definite matrix thus it can be decomposed as follows 

BB T = V diagip 2 p 2 ,p 2 } v T 

2 2 2 2 2 2 
where diag{PpP 2 ,P 3 ) is a diagonal matrix whose elements are ppp 2 and P^. 

Consequently 


L = ±(BB T ) 1/2 = V diag{±p l5 ±P 2 , ±P 3 ) V 1 


Substitution of (12) into (11) yields 


A q = T(BB T )' 1/2 B = V diag{±Pj, ±P 2 > ±P 3 ) V T B 


( 12 ) 


(13) 


It can be verified that to obtain maximum of p*(A) we need to choose the plus signs 

1 1/2 

in (13). We designate it by choosing the plus sign in front of (BB ) ; that is 


(14) 



which is the sought solution of Wahba’s problem. 

The expression given in (14) is also the solution of another problem as 
discussed next. 


III. THE CLOSEST ORTHOGONAL MATRIX 

Consider the following problem. Given a real matrix, B, what is the closest (in 
the Euclidean-norm sense) orthogonal matrix to it? To solve this problem denote the 
square of the Euclidean norm of the difference between B and any same order real 
matrix, A, by s(A); that is 

s(A) = ||B - A|| 2 

(where | | . | | denotes the Euclidean-norm) and find the 3x3 orthogonal matrix, A, which 
minimizes s(A). It can be easily shown that 

s(A) = tr[(B - A)(B - A) T ] 


thus 


s(A) = tr(BB T - BA T - AB T + AA T ) 

Using the fact that A has to be orthogonal and the properties of the trace operation 
it can be easily shown that 


s(A) = tr(BB T ) + 3 - 2tr(AB T ) 


(15) 


Obviously, that A which minimizes s(A) is the A which maximizes the term tr(AB 1 ). An 
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inspection of (2) reveals that this particular A is also the solution to Wahba’s 
problem. This result can be stated as follows. The closest orthogonal matrix to B, 
where B is as defined in (3), is the solution to Wahba’s problem. Indeed if we 
proceed with finding that orthogonal A which minimizes (15), we will obtain the 
result given in (14); namely, 

A q = (BB T )' 1/2 B (16) 

Consequently, any solution to the closest orthogonal matrix problem is also a 
solution to Wahba’s problem. This conclusion will be exploited in the ensuing. 

The solution expressed in (16) to the closest orthogonal matrix problem was 
obtained and investigated quite extensively in the past [12 - 19]. The solution of A q 

using (16) is cumbersome. Various iterative solutions have been investigated [15 - 

i9]. 

Another solution to the closest orthogonal matrix problem, and hence to Wahba’s 
problem, makes use of the singular value decomposition (SVD) of A q . This solution is 

presented next. It is well known [20] that any matrix, and therefore also B, can be 
decomposed as follows 

B = USV T 

where U and V are 3x3 orthogonal matrices and S is a diagonal matrix whose elements 

T 

are the nonnegative square roots of the eigenvalues of B B. It can be shown that 

T 

A = UV 1 
o 

The latter was used in [21] to solve Wahba’s problem. 

IV. POLAR DECOMPOSITION 

It is well known [22] that B can be decomposed as follows 


B = PH (17) 

where P is orthogonal and H is symmetric. This decomposition is known as polar 
decomposition (PD). It was shown [23] that P is precisely the orthogonal matrix 
closest to B; that is, P of the polar decomposition is the solution to Wahba’s 
problem when B is as defined in (3). We can write therefore 

B = A H 
o 

(where A q is, of course, the closest orthogonal matrix to B). This yields 

A = BH' 1 (18) 

o 

We wish now to utilize the PD concept for solving Wahba’s problem. We consider two 
cases as follows. 

IV. 1: The Error Free Case 

Assume now that both sequences of vectors and r- i=l,2,3,...K are error free. 
We can then write b-=Ar. Substitution of this equation into (3) yields 
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(19) 


K 


K 


K 


B = V abrf = Y a-Ar.r* = A y a-r.r: 
L 111 L 111 L 111 

i = 1 i = 1 i = 1 

Define now the matrix R as follows 

K 


R = y a.r.r' 
L ill 

i = 1 


then from (19) we obtain 


B = AR 


( 20 ) 


( 21 ) 


where R is a symmetric matrix. Comparing (21) with (17) it is easy to see that in 
this case (21) is the PD of B where A=P and R=H. It is clear then that A q =A. In this 

case A q can be found as follows 


A 0 = BR' 1 (22) 

provided that in constructing R, according to (20), we use at least 3 non-collinear 
vectors. (This assures that R is invertible.) 

IV. 2: The Actual Case 


In practice the vectors bj are contaminated by measurement noise. However, since 

the position of the body and the time of measurement are known within a high degree 
of precision, the error in the determination of the r. vectors is negligible. Denote 

the error in bj by n. then we can write that 

b. = n. + Ar- 
i i i 

Using the last equation in (3) we obtain 


This can be written as 


B = 


l a i ( "i + Ar i )r i 
i = l 


which yields 


K T K T 

B = y a-n-r- + y a-Ar-r: 

L ill L ill 

i=l i= 1 


K T K T 

B - y anr: = A Y a-r.r. 

L ill L ill 


i = l 


i = l 


Using (20) we obtain from the last equation 
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(23) 


A = 


K 


B - l 

i = l 


T 

anr • 

1 1 i 



We can now use the last equation to obtain the "best" estimate of A. We note that B 

A 

contains all measured information, therefore we compute A, the "best" estimate of A, 
as the conditional expectation of A given B [24]. Performing the conditional 
expectation on both sides of (23) yields 


E(A/B) = 


K 

IB- I 

i = 1 


ajElm/Bli-y 


R 


-1 


It is assumed that the measurement errors are unbiased, therefore 


(24) 


Efm/B) = 0 (25) 

(The latter assumption is based on the premise that the measurement biases have been 
removed or else are very small. If this is not the case, there is no way to obtain 
the correct attitude from the biased measured vector no matter what algorithm is 
used.) Substitution of (25) into (24) yields 

E(A/B) = BR' 1 


thus 


A 

A = 


BR 


-1 


(26) 


where B and R are computed according to (3) and (20) respectively. 

Note that this result was first obtained by Brock [13, eq. (5)] in a way 
unrelated to the notion of polar decomposition and with no consideration of the 
randomness of n. 

If n- are very small or the number of measurements is large such that the 
particular realization of n. has a negligible mean, which complies with the 

1 A 

assumption in (25), then the computation of A according to (26) yields an accurate 
estimate of A. When this is not the case, the estimate can be quite erroneous. It is 

interesting to note that when K<4, A zeros the cost function of Wahba’s problem which 
is given in (1) as follows 

p(A) = j l aj |bj - Arjl I 2 
i = 1 


even if A is not equal to A. This is a result of the approximation b-=Ar- which was 

A A , 1 . 1 A 

made in the derivation of A. However, while A drives p(A) to its minimal value, A is 
not necessarily orthogonal. (Recall that we seek the orthogonal matrix which 

A 

minimizes p(A)). We can correct the non-orthogonality of A by the application of one 
orthogonalization iteration as follows [17, 18] 

A’ = 0.5(A" T + A) (27) 

A 

This operation yields a close to orthogonal matrix. A’, which is usually also closer 
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to A. (The superscript -T denotes the inverse of the transpose.) We can, of course, 
bypass the computation of A by using (26) in (27) to obtain 


A’ = 0.5(B’ T R + BR' 1 ) 


(28) 


V. NUMERICAL COMPARISON 

Five possible solutions to Wahba’s problem are considered as follows. 

(1) QUEST 

Use the algorithm QUEST [6] to obtain, q, the quaternion which corresponds to 
the solution matrix of Wahba’s problem, and then use q to compute the solution 
matrix itself which we denote by A 

(2) ITERATIVE ALGORITHM (IA) 

Apply the iterative orthogonalization algorithm [17, 18] starting with the 
computation of B according to (3) and then continue with 

A q = B (29.a) 

A A_T A , v 

A j+1 = 0.5(A j 1 + A.) (29.b) 

which converges to the solution of Wahba’s problem given in (13). We denote the 
final matrix by A^. 


(3) SINGULAR VALUE DECOMPOSITION (SVD) 
Apply the SVD algorithm to decompose B into 

B = USV T 


and compute 


^svd 


T 

= UV 


(30.a) 


(30.b) 


As explained in Section III, A gv( j too is the solution of Wahba’s problem. 


(4) FAST OPTIMAL MATRIX ALGORITHM (FOAM) 

Use the FOAM algorithm [25] to obtain the solution matrix to Wahba’s problem. We 
denote the computed solution by A^ om . 


(5) POLAR DECOMPOSITION (PD) 

Compute the matrices B and R, the latter according to (20), and then calculate 
the estimate of the solution to Wahba’s problem according to (26) 
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(26) 


(6) IMPROVED POLAR DECOMPOSITION (IPD) 

Compute an improved estimate of the solution to Wahba’s problem by performing 
one orthogonalization iteration on the preceding estimate. The overall algorithm 
is as in (28) 

A’ = 0.5(B’ T R + BR' 1 ) (28) 


V.l Results with Simulated Data 

The five algorithms were tested with simulated data. The importance of tests 
with simulated data stems from the fact that using real data we do not know the 
correct attitude. This constitutes a major difficulty since the difference between 
algorithms may be smaller than the difference between the correct attitude and the 
computed ones. Only when we use simulated data can we observe the difference between 
the computed attitude and the correct one. The simulated measurements of vectors in 
body axes were obtained by transforming the reference, i*j, vectors to body axes using 

A, the correct attitude matrix, addition of a noise component to each component of 
the transformed vector and normalization of the resultant vectors. The added noise 
components had a zero mean and a standard deviation value of 0.144. Typical 
simulation results are shown next for four and three measured vectors. Three cost 
values were computed in order to evaluate the accuracy of the results. The cost p is 
Wahba’s cost function computed according to (3) for the particular solution matrix. 
The cost f is the Euclidean norm of the difference between the particular solution 
matrix and the correct attitude matrix. Finally, the cost J is a measure of the 
non-orthogonality of the solution matrix. It is the Euclidean norm of the matrix 
T 

XX -I where X is the particular solution matrix. 

V.1.1 Four reference vectors 

'.2672611 [-.6666671 [ .2672611 [-.447214' 

r, = .534522 r~ = -.666667 r~ = -.801784 r . = .894427 

1 [.801784J 1 L-.333333J [ -534522J [ .000000 

Four "measured" body vectors 

'.8153991 [-.8722141 [.2902031 [-.118959' 

b, = .577901 b 0 = -.075280 b^ = -.206009 b. = .679197 

1 -.033975 .483296J J [ -934528] L'- 706940 


Four weights 

aj = .100000 a 2 = .300000 a 3 = .400000 a 4 = .200000 

The correct attitude matrix 

• .764744 .293558 .573576' 

A = -.636031 .486370 .599090 

-.103103 -.822963 .558660 
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Solutions 


A = 
qst 




A 

fom 


A 

A = 


A 

A’ = 


‘ .770135 

.274589 

.575754' 

p = .44078E-03 
*qst 

-.629265 

.474889 

.615228 

f t = .37578E-01 
qst 

-.104485 

-.836111 

.538518 

J , = .20588E-06 
qst 

‘ .761290 

.266299 

.591204' 

p- = .22933E-03 
Mtr 

-.639697 

.457436 

.617689 

f. f = .67264E-01 
itr 

-.105948 

-.848432 

.518593 

J. t = .19037E-06 
itr 

' .761290 

.266299 

.591204 

p syd = .22933E-03 

-.639697 

.457436 

.617689 

f , = .67264E-01 
svd 

-.105948 

-.848432 

.518593 

J , = .18014E-15 
J svd 

’ .770135 

.274589 

.575754 

p fom = • 44078E -° 3 

-.629265 

.474889 

.615228 

f, = .37578E-01 
fom 

-.104485 

-.836111 

.538518 

J r = .86667E-16 
fom 

' .770556 

.263174 

.561689' 

p = .67846E-04 

-.654729 

.455528 

.628148 

f = .75595E-01 

-.143061 

-.851596 

.551271 

J = .11190E+00 

' .768038 

.263582 

.584080' 

p’ = .33350E-03 

-.630931 

.473739 

.615274 

f’ = .52240E-01 

-.115625 

-.840565 

.530461 

J’ = .25319E-02 


V.1.2 Three reference vectors 


.2672611 


[-.6666671 


.2672611 

.534522 

r 0 = 

-.666667 

^ 'J ~~ 

-.801784 

.801784 

z 

-.333333 

J 

.534522 


Three "measured" body vectors 


.8153991 


[-.8722141 


.2902031 

.577901 

b 9 = 

-.075280 

b 3 = 

-.206009 

-.033975 

z 

.483296 

J 

.934528 


Three weights 

a, = .125000 a 2 = .375000 a 3 = .500000 
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The correct attitude matrix 


A = 


’ .764744 
-.636031 
-.103103 


.293558 

.486370 

-.822963 


.573576' 

.599090 

.558660 


Solutions 





' .767731 

.276451 

.578069' 

Pqst 

f t = 
qst 

^qst ~ 

.55884E-03 

-.631488 

.479442 

.609392 

.29705E-01 

-.108683 

-.832893 

.54265 8_ 

.67617E-07 

‘ .758264 

.271018 

.592946' 

Pitr = 

.23600E-03 

-.643834 

.454336 

.615676 

f itr = 

.67219E-01 

-.102537 

-.848604 

.518997 

J itr = 

.32845E-06 

’ .758264 

.271018 

.592946' 

Psvd - 

: .23600E-03 

-.643834 

.454336 

.615676 

f svd = 

: .67219E-01 

-.102537 

-.848604 

.518997 

^svd " 

: . 1 1 102E-15 


^fom 


' .767731 
-.631488 
-.108683 


.276451 

.479442 

-.832893 


.578069' 

.609392 

.542658 


p fom = - 55884E -° 3 
f fom = • 29705E -° 1 
J fom = -30626E-15 


' .739265 

.275664 

.586784' 

p = .11783E- 14 

-.664499 

.459428 

.635984 

f = .97131E-01 

-.172692 

-.839769 

.575035_ 

J = .16640E+00 

' .753716 

.268839 

.600058' 

p’ = .60457E-03 

-.645610 

.483007 

.593789 

f’ = .53687E-01 

-.131702 

-.833708 

.539069 

J’ = .53638E-02 


We observe that, as expected, Aj and A sv( j are practically identical. We also 
observe that as expected, for three measured vectors (K=3) Wahba’s cost, p, for A is 

A 

practically zero. The single normalization cycle which generates A’ improves the 
orthogonality (reduces J) considerably. This A comes at the expense of an increase in 
p. For four measured vectors (K=4), p for A is similar in value to that of the other 
algorithms, and again, the single normalization cycle improves orthogonality 
considerably at the expense of Wahba’s cost. 
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V.2 Results with UARS Data 

The following are results of the application of the five algorithms to data 
measured on-board the Upper Atmosphere Research Satellite (UARS). UARS was deployed 
on September 15, 1991 at 04:23 GMT by the shuttle spacecraft Discovery which was 
launched on September 12, 1991 , at 23: 12 GMT. The data were measured on September 30, 
1991 at 18:32:31.206749916 GMT. The first vector corresponds to the Sun Sensor, the 
second to the triad of Magnetometers, and the third to the Infra-Red Horizon Sensor. 

The reference vectors 

'-.992324] [-.814177] [ .543295' 

r, = -.113458 r 0 = .550862 r~ = -.542620 

1 [-.049192J z [-.183487J [.640619 

The measured body vectors 

'-.810765] [-.455867] [ .002528' 

b, = -.294952 b~ = .186491 b^ = .003031 

[-.411 403 J [-.870291J [ .999992 

Three weights 

aj = .243291 a 2 = .002506 a 3 = .754203 

Solutions 

’ .826549 .178850 -.5336941 p q$t = .96423E-03 

A = .178119 .816336 .549426 

qst 

.533939 -.549189 .642885 J , = .20183E-06 

L J qst 

‘ .832537 .172669 -.526372] pj = .89246E-03 

Aj = .180280 .814010 .552166 

[ .523814 -.554593 .646564] Jj = .15599E-07 

' .832537 .172669 -.526372] P svd = .89246E-03 

A syd = .180280 .814010 .552166 

.523814 -.554593 .646564 J . = .19700E-15 

L J svd 

' .826549 .178850 -.533694] p fom = .96423E-03 

A fom = - 178119 -816336 .549426 

[.533939 -.549189 .642885] J fom = .46516E-15 
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‘ .818163 

.211577 

-.510709' 

P = 

.1361 IE- 12 

.182985 

.778550 

.508996 



.466235 

-.700019 

.572641 

J = 

.28575E+00 

‘ .837978 

.176650 

.517931' 

P’ = 

.39206E-02 

.217866 

.767101 

.613477 



.503300 

-.622646 

.606208 

J’ = 

.16305E-01 


Here too we observe the identity between and A sv( j. As before, we also 
observe the reduction in J at the expense of an increase in p when a single 

A A 

orthogonalization cycle is applied to A to generate A’. 

V.3 Time Consumption Analysis 

A computation-time measurement was performed on all five algorithms using the 
simulated three and four measured vectors. The runs were made on a VAX 9210 computer 
employing the VMS Version 5.4-2 operating system. The time measurement routine used 
the internal machine clock at a resolution of 10 msec. In order to increase the 
resolution, the runs were performed over 50000 successive solutions and the total 
time was then divided by 50000. The results are presented in Table I. 

Table I: Algorithm Computation Time (msec). 



A , 

qst 

A. 
i tr 

^svd 

^fom 

A 

A 

A 

A’ 

Three 

measured 

vectors 

0.0890 

0.790 

0.548 

0.060 

0.058 

0.084 

Four 

measured 

vectors 

0. 1060 

0 . 694 

0 . 526 

0.070 

0.068 

0.094 


Note the decrease in computation time of A^ when the number of measured vectors 

increased from 3 to 4. This is due to the fact that in the four vector case the 
convergence criterion was met after only 7 iterations whereas in the 3 vector case 8 
iterations were performed until the same convergence criterion was met. In all our 
tests it was found that when a fourth measured vector was added, less iterations were 
required. This stemmed from the fact that when a fourth measured vector is added the 
orthogonality of B increases provided the fourth vector is not a linear combination 
of the other three. The decrease of the computation time of A gv( j with the increase of 

the number of measured vectors is not consistent. It depends on the number of 
iterations needed for the completion of the SVD calculations. 
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VI. CONCLUSIONS 


It was shown that the solution to Wahba’s problem is the closest orthogonal 
matrix to B where B is defined on the measured vectors and on weights associated with 
their measurements. The weights signify the confidence assigned to the measurements. 
The matrix B includes all the information contained in the measurement. 

Once it was established that the sought solution is the orthogonal matrix 
closest to B, algorithms for computing that orthogonal matrix were considered, and an 
algorithm was discussed which is based on the polar decomposition of matrices into 
the closest unitary (in our case: orthogonal) matrix and a Hermitian (in our case 
symmetric) matrix. The accuracy of the algorithm increases with the accuracy of the 
measurements and with their unbiasedness. If the measurements are error free the 
algorithm yields the exact solution. 

When only three measured vectors are used the new algorithm yields a solution 
which zeros Wahba’s cost; however, the solution is not necessarily orthogonal. An 
application of one orthogonalization iteration to the solution matrix constitutes a 
modified algorithm which yields a better solution. Although the latter algorithm 

generates a matrix which increases Wahba’s cost. The new matrix is closer to 

orthogonality. We note that the same iteration cycle if applied repeatedly to B 
itself, yields eventually the optimal solution as shown in the examples; however, 
since B is usually quite far from orthogonality, it takes several iterations to 

obtain the solution. 

The advantage of the algorithm is in its simplicity which enables its use for 
obtaining first cut solutions using "back-of the envelop" like programs such as 
MathCAD. Another advantage of the first new algorithm is its ability to indicate the 
precision of the measurements. This stems from the fact that generally the closeness 
of the solution matrix to orthogonality is indicative of the precision of the 
measurements. It is interesting to note that the fact that the two PD algorithms 
yield the exact solution in the noise-free case is analogous to the fact that the 
largest eigenvalue of the 4x4 K matrix used in the QUEST algorithm is precisely 1 in 
the noise-free case. 

The two PD algorithms were tested vis-a-vis other popular algorithms using 

simulated and real UARS data. 
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